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The position B vg of a field-tuned superconductor-insulator quantum transition occuring in 
disordered thin films is examined within the mean field approximation and starting from a 
hamiltonian of BCS type (favoring the s-wave pairing) with a random potential, and effects of 
an electron-electron repulsion on the transition field B vg are also considered. Our calculation 
shows that the (microscopic) disorder-induced reduction of B vg suggested experimentally cannot 
be explained without taking account of the familiar interplay between the randomness and the 
electron-electron repulsion enhancing the quantum superconducting fluctuation in such systems. 
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§1. Introduction 

There have been many reportsBSi) on resistive data suggestive of a strong quantum fluctuation 
effect in disordered thin films under magnetic fields perpendicular to the film plane. In such systems, 
the temperature dependence of the sheet resistance R(T) in higher fields than a field B c is typically 
insulating, while in lower fields it decreases on cooling, in most cases, in a manner suggesting the 
presence of a magnetic field-tuned superconducting transition at zero temperature (T = 0). It 
was found experimentallyHi) that B c decreases with decreasing the film thickness or increasing the 
impurity density, both of which enhance the high temperature sheet resistance R r . 

Fisherft has proposed by extending the theoretical argument for zero field case& to nonzero 
fields (B y£ 0) that, at a field B vg implying such a quantum superconducting transition point 
to be identified with B c , the resistance curve R(T) at lower temperatures should be independent 
of T and take a universal constant value. However, extending directly the description!) in zero 
field to B case is theoretically invalid because a superconducting phase in the presence of 
the field-induced vortices is not the Meissner phase but a kind of a vortex glass (VG) phasJli^ 
induced by a pinning disorder and with no Meissner effect. On the other hand, one of the present 
authors has pointed outH' that, in the field range where the pinning disorder is negligible, the 
insulating behavior of resistance in T ^ can intrinsically arise due to the quantum superconducting 
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fluctuation, and that, if the normal contribution a n to the conductivity is negligible, a flat R(T) 
curve with approximately the value of quantum resistance Rq = 7rh/2e 2 = 6.5(kJ7) is expected 
along the quantum melting (crossover) line at low temperatures which, as far as the ordinary 
dirty limit is valid, is insensitive to T. Also, it has been recently!) pointed out in terms of a 
phenomenological Ginzburg-Landau (GL) action with pinning disorder terms that, irrespective of 
the details of dynamics of the superconducting fluctuation, the disorder-induced contribution to 
the T = conductance at B vg should be, in contrast to the original argument,!) a nonuniversal 
constant dependent on the strengths of pinning and fluctuation. 

However, in the present problem including a static disorder on the electronic level, a simple 
analysis starting from a GL action is insufficient, and more or less, one needs to return to a 
microscopic electronic hamiltonian. In a previous work,!) we have found, by deriving a quantum 
GL action, a significant enhancement of quantum fluctuation in nonzero fields due to an interplay 
between the microscopic disorder and a repulsive interaction between electrons. However, the 
familiar derivation of a GL action done there has neglected spatial variations of the coefficient of 
each term in the resulting GL model. A model of such spatial variations used widely in a GL 
theory is a random potential for the pair- field implyingBS' a randomness of T c o, where T c q is the 
mean field transition temperature in zero field and in clean limit. To clarify whether the ideal!) 
based on the quantum superconducting fluctuation is applicable to real systems, we have to take 
account of such a GL random potential term which induces, in a finite field, a pinning disorder0) 
for the vortices. A crucial point is that, in the present issue, such a GL random potential arises 
from the microscopic disorder u(r) which simultaneously enhances the (quantum) superconducting 
fluctuation. Note that, for a decrease and vanishing of resistivity in B ^ 0, the superconducting 
fluctuation effect and the vortex pinning are competitive with each other. Due to this competition 
originating from the same microscopic disorder, the disorder-induced decrease^!) of B vg may not 
be necessarily expected theoretically. The purpose of the present study is to give an answer to this 
question by starting from a microscopic model. 

To examine the quantum superconducting fluctuation!) and the vortex glass fluctuation,!) arising 
from the pinning disorder, on the same footing, one needs an appropriate quantum GL action with 
spatially-dependent coefficients. Within the approximation taking account only of the pair-field 
(superconducting order parameter) \& belonging to the lowest Landau level (LLL), the quantum 
GL action is phenomenologically expected to take the following form 



S ran = J d 2 r /3]>>(0)+7M)|*Mr)| 2 + J du ( fy(r)|tf(r, u)\ 2 + ^|tt(r, u 



(1.1) 



where u is an imaginary time, u a Matsubara frequency, 5fj,(r) a static random potential implying 
a spatial variation of T c q, ty(r,u) (= A rl//2 (0) < ip-\(r, u)ip i(r, u) >) = J2w ^u( r ) exp(— iwu), and 7 
measures!^ the dissipation strength and was denoted as 71 in ref.9. Since we focused on the LLL 
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modes, eigenvalues of gradient terms were already absorbed into //(0). For simplicity, a randomness 
of the quadratic term was taken in eq.(l.l) in a local form, i.e., up to the lowest order in the gradient. 
In §2, an action corresponding to a replicated version of the action (1.1) will be microscopically 
derived by neglecting the electron repulsion and focusing mainly on the low T and high B region. In 
§3, our formulation for determining the Gaussian vortex-glass transition field -B V g,o at T = will be 
explained and is applied both to the ordinary dirty limit with no repulsive interaction and to the case 
with an electron-repulsion. We find that, in contrast to the experimental observation ,§'§!) -B V g,o in 
the ordinary dirty limit increases with increasing the (microscopic) disorder, mainly reflecting an 
enhancement of the corresponding H C 2(0) (= H^ 2 (0) oc T c q/t). Further, the analysis is extended to 
the case with a repulsive interaction, and we find that the quantum LLL fluctuation and a decrease 
of H C 2(0) due to an inclusion of an electron-repulsion are origins of a reduction of B vg consistent 
with experimental observation. An attempt of computation on low T behaviors of 7 is given in 
§4 to demonstrate the presence of a region in which the GL coefficients are insensitive to T. Our 
results are summarized in §5 together with a comment. 

§2. Model and Derivation of Pinning Vertex 

As a microscopic basis of our analysis, we first start from the s-wave BCS Hamiltonian with 
random potential in two dimensional (2D) case 



H 



£4(r) I ^-HV - ^A(r)) 2 + «(r) ) ^(^+^^(1)^(^(^(1 



2m 



(2.1) 



For the moment, we will not take account of electron-electron repulsions and postpone its inclusion 
to §3. In eq.(2.1), gees < 0, </>o is the flux quantum for Cooper pairs, ip a is an annihilation 
operator of electron with spin a (=], |), and the applied field B = curlA is perpendicular to the 



2 D plane. The random potential u(r) has zero mean and the Gaussian ensemble u(r)u(r') = 
(2ttN(0)t)~ 1 5 ( - 2 \t — r'), where r is the elastic relaxation time, iV(0) the density of states at Fermi 
surface, and the overbar denotes the microscopic random average. Following the study of zero field 
transition point in the literature,© the dependence of physical quantities on the film thickness is 
assumed to arise entirely through the repulsive interaction effect. In most part of this paper, our 
analysis is based on the ordinary dirty limit in 2D (with no interaction effect) in which the disorder 
strength is measured by {kpl)~ l = t/(2itN(0)1 2 ). 

In our microscopic analysis which follows, the so-called quasi-classical approximation will be used. 
In this treatment, a gauge-invariant gradient Q = — iV + 2ttA/(J)q operating on a f is regarded as 
a wavevector (c-number) in the electronic process. Further, before focusing on the LLL modes of 
\E', the coefficients of the quadratic and quartic terms in the action are found to be functions of 
Q 2 , while Q 2 is transformed, after operating on a LLL ^ mode, into r^ 2 = 2ttB/<Po where B is the 
uniform flux density. Hence, after this transformation Q 2 — ► r^ 2 , the coefficient of each term in the 
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action, except a term leading to the vortex pinning, may be written as a constant. In addtion, we 
note regarding the validity of the LLL approximation that, near T = 0, the field range over which 
the LLL approximation for the pair field is valid is quite broad even below H C 2(0) according to 
ref.8. Under these conditions and within the ordinary dirty limit neglecting the localization effect 
and the electron-electron repulsion, the coefficients n(0), 7, and C/4 in the resulting eq.(l.l) are 
well-known and given byi© 



U, 



"-mug™* 



M C°)(0) = (iV(0)| 5B cs|) _1 - 47rr/3- 1 Er(2e;S) 



e>0 



,(0) 



47T/r 1 r 2 E(r(2e;£) 



(2.2) 
(2.3) 
(2.4) 



e>0 



where T(2e; B) = (2|e| + 2-kBD / '(^o)' 1 ^ 1 implies r Q (2e) with Q 2 = 27r£/</> , T q (e + e') = r _1 (|€ + 
e'| + Dq 2 )^ 1 is the diffusion propagator with momentum q, e and e' are Matsubara frequencies 
for fermions, and the diffusion constant D is given in terms of the mean free path I = kpr jra by 
D = 1 2 /(2t). Through this paper, the (impurity-averaged) electron propagator is, as usual, given 
by G k (e) = (ie - £ k + (i/2r)sgn(e))- 1 . 

Next, to explain how the Sfj,(r) term, neglected in previous studies,^ appears in the present 
formulation, let us first imagine the phenomenological action (1.1) to be rewritten in a replicated 
form. Following the standard treatment for an action of GL type, the replicated form© of eq.(l.l) 
will be given by 



5r, 



E 

a=l 



/3E(M0) + 7M)l*L Q) (r)| 2 + y / du Jj P i a \u)\ 



tJ dUl I dU2 i 



a'=l 



U p (k) ( a ) („/) 



(2.5) 



where a and a' denote the replica indices, U\ and u 2 are imaginary times, and the bare (It- 
dependent) pinning vertex function U p (k.) and p^\u) are Fourier-transformations, respectively, 
of £/i(r) fyi(r') and of (^(r, u)\ 2 . As seen below, the k-dependence of the bare pinning vertex 
U p (k) is characteristic of the high B and low T region. 

A microscopic derivation of the last term of eq.(2.5) will be explained here. Within the quasiclas- 
sical approximation ,0) the third (replica off-diagonal) term of eq.(2.5) generally takes the following 
form 

As? an = -yE E / ^EE-niQ*]; <^0 

a=l a'=l w ui' 



x(*(f)( ri )r ^\v 2 ) (^:\v,)f *^(p 4 ; 



(2.6) 



ri,r2,r 3 ,r 4 - 
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where Q s = — iV s + (27r/^ )A(r s ), JF([Q S ]) is a bare vertex of the quartic term induced by the 
random-average, and [Q s ] implies dependences on QJ, Q2, Q3, and Q4. Eq.(2.6) includes terms 
off-diagonal in the replica indices as a reflection of the fact that no impurity lines carry finite 
frequencies. Due to the use of the quasiclassical approximation, the spatial nonlocality, arising from 
the nonzero magnetic field, in T(- ■ •) appears only in the form of the gauge-invariant gradients Q s 
operating the pair-fields *(r s ). Hereafter, frequency dependences of T will be neglected. 

Below, let us see how AS 1 ™^ will lead to the form of the last term of eq.(2.5) in a consistent 
approximation with the derivations of ll[ \ /U^(0), and 7^. Namely, as in the derivation of 
Uf \ the dia grams expressing T are selected under the condition, equivalent to the neglect of 
localization effect of noninteracting electrons, that the impurity lines, each of which carries the 
factor (2ttN(0)t)~ 1 , do not cross to each other in the diagrams. The diagrams thus obtained 
expressing T are the same as those in ref.14, where the diagrams were examined from the viewpoint 
of mesoscopic fluctuation, and their examples belonging to the same family are given in Fig.l. 
Effects of an electron-repulsion on AS™ an will not be examined here. We simply note that, just like 
those in U4, each term perturbative in a short-ranged repulsive interaction in the function JF([Q S ]) 
is convergent in T — > limit. The first three diagrams of those in Fig.l are given in Fig. 2, on which 
we will focus here to clarify key features common to all diagrams in Fig.l. In the figures, the single 
solid curves denote G&(e)'s, four double solid lines denote ^(r.;) or **(rj), and Qj (Q*) operates 
on *(r.j) (\I/*(rj)). The sum of diagrams in Fig.2 contributes to JF([Q S ]) in the manner 

^ 2 >([Q a ]) = (tf(o)/r 2 ££r QI (2^ r q+Q ,_ Q2 ( e + e ') 

€>0e'>0 q 

xr q ( e + ^ G k+Q2 ( e )G k (-e)G k+Q * +q (-e')G k+q (e') (1^ + 1^ + 1^), (2.7) 

where 

I (2a) = / k G k+QI ( e )G k (- e )G k+q+Q4 (- e 0G k+q (e'), 

/(») = -(27rriV(0))- 1 [ G kl+Q j( e )G kl (- e )G kl+q (e') / G k2+QJ ( e )G k2+q+Q4 (-e')G k2+q (e / ), 
and 

= -(27rriV(0))- 1 / G kl+Q j( e )G kl (- e )G kl+q+Q4 (-e') / G k2 (-e)G k2+q+Q4 (-e')G k2+q ( e / ), 

Jki Jk2 

(2.8) 

and /( 2a ) is the contribution of Fig.2(a) and so forth. However, a cancellation occurs between the 
three diagrams in the manner that the sum 

/(2a) + j(26) + 7 (2) reduces to 27TiV(0)T 4 (2e + 2e' + 
D(Ql) 2 + DQl + 2D(q 2 + q . (Q 4 - Q*))). It is easily understood that one more cancellation similar 
to this arises when all diagrams of Fig.l are summed up, and the contribution of Fig.l to T becomes 

^)(Q,AQ)= ^^jj EE( r Q( 2f ) r Q( 2£ ')) 2 ( r ^( £ + Or q (e + 
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x (e + e' + DQ 2 + Dq 2 + Dq • AQ) V, (2.9) 

where AQ = Q4 — Q\ = Q3 — Q2, and all Q 2 and (QJ) 2 were expressed as Q 2 , because the 
external pair-fields are assumed to be in LLL. Further, using the fact that, after operating a 
LLL eigenfunction, Q 2 changes into r^ 2 = 27rB/<f)o, eq.(2.9) implies that J 7 ^ takes the form 
•F«(AQ) = U p fW(t; r B AQ), where 

\ 2 

4 rgr 



and t = 8ttt/ (/3l 2 r B 2 ). Although other families of diagrams are also found to have similar structures 
to JFW, the full expression T obtained after summing them up has a highly complicated AQ- 
dependence. Here we merely mention that, if formally setting AQ = 0, the full T becomes the 
simplified form: 

j-(aq = o) ~ 12 j £ £ r Q (2 e )r Q ( 2e ') jT ( r q ( e + e ')) 2 , (2.11) 

which is estimated as 0.14U P in low t limit. Since, as is seen below, it is difficult to write down 
a concrete expression of U p (k) resulting from the full JF(AQ), we will merely explain below some 
properties of the full f7 p (k). 

At this stage, let us rewrite the terms of eq.(2.5) entirely in terms of the LLL fluctuation field 
ipo^ip, ui) so that p£\u) = J2w P^(^) w)e -lwu is expressed in a Landau gauge, by 

p^(k, fi) = ^ exp(-^i + ipk x r 2 B ) £ + V2, " + fi))* 4% - ky/2, w), (2.12) 

p ' a; 

where ip^ip, uj) is defined as ^(r) = ^ p <A)(Pj ^)^o,p( r ) in terms of a LLL eigen function uo lP - If 
the AQ-dependences in T are neglected, the resulting pining vertex is local (i.e., the function U p (k) 
is k-independent). This approximation is valid in high T and low B region defined by t 3> 1, i.e., 
T S> 0.15T c o5/i^c2(0). In contrast, at low T and high B of our interest, any microscopic length on 
the pair-fields is scaled by re, and thus, the AQ-dependences are no longer negligible. Actually, 
we have the relations in a Landau gauge 

^AQ(M/M( ri ))^ (Q) (r 2 )(^( a ')(r 3 ))*vI/(«') (r4) | ri)r2)r3)r4 _ 

i W(z x k) ^-*>')^| ( ^) ( p _ ky/2)y^{p + V2)(rf*V + *y/2)) Vj*V " *y/2), 
( A Q) 2 (^W( ri ))*vI/( a )(r 2 )(*( Q ')(r3))^(«') (r4 )| rii 



p,p' 



,r2,r3,r4^r 



(2.13) 



p,p 



6 



and so on, where t> q = exp(— q 2 r^/2). Namely, any rgAQ-dependence changes, after the Landau 
level representation, into a r#k-dependence or a constant. Hence, U p (it) in eq.(2.5) generally has 
the form 

U p (k) = U p foo(t; kr B ). (2.14) 

If the Fermi surface is isotropic, /oo is a function of k 2 r 2 B . Below, for brevity, we will focus on this 
isotropic case. At low temperatures t <C 1, /oo is not sensitive to the material parameters included 
in the definition of t but depends merely on k 2 r\. Unfortunately, it is difficult to find an exact form 
of such /c-dependences in terms of eq.(2.13) which implies that a AQ-expansion does not reduce to 
a systematic k-expansion. Nevertheless, we have tried to estimate the ^-dependence of /oo at T = 
occuring from the low frequency limit of frequency- integrals, because the corresponding dependence 
occuring from higher frequencies is clearly regular, and a large /^-contribution is cut off by the i>k- 
factor appearing everywhere in the LLL diagrammatics. Using the relations (2.13), a nonanalytic 
/c-dependence arising from the lowest frequencies is found to be less singular than krs ln(/crs) _1 
which vanishes in k — > 0. We note that, since k is always accompanied by rg, the pinning strength 
occuring after a k-integral in a physical quantity is measured by U p /rg ~ (Eft)~ 2 in the dirty 
limit in t < 1. Since the quantum fluctuation strength U^/{r 2 3 r y) becomes, as seen below, of the 
order (Epr) in the dirty limit, this implies that the relative (effective) pinning strength Upj/U^ 
is 0((Eft) _1 ). It is verified in terms of eqs.(2.2) and (2.11) in t 3> 1 limit that this is of the same 
order as the thermal counterpart, denoted as A e fj in ref.15, (3U p (t S> l;k = 0)/U4 in the dirty 
limit. Therefore, in the ordinary dirty limit, the vortex pinning strength relatively increases at any 
T with increasing disorder with strength (Epr)^ 1 . 

§3. Quantum VG Transition Point 

In this section, we examine the T = VG transition point -B V g,o i n the Gaussian (i.e., mean field) 
approximation using the LLL-GL action in the dirty limit derived in §2, and, based on this result, 
effects of the electron-electron repulsion on -B V g,o are qualitatively studied. The VG transition 
field is dennedBHEl) by examining where the uniform VG susceptibility Xvg (defined below) at zero 
frequency tends to diverge. For a moment (until reaching eq.(3.16)), the GL action is assumed to 
be well-defined in low T limit. 

The key quantity Xvg in this section is proportional to / f \ < 1 I' w= o(ri)^'* =0 (r2) > | 2 and, 
consistently with eqs.(2.5) and (2.11), is given by 

xv g = £ W&KW, (3.i) 
p-p' 

where G(p,p', u>) =< (po(p, u>) ty?o(p'> w ) > ^ s the LLL fluctuation propagator defined prior to the 
random averaging. On the other hand, the renormalized 2D LLL superconducting fluctuation, 
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defined through the random-averaged (replica-diagonal) fluctuation propagator 

g(u) = g(p,p>, lo) = ( mo) + tM + s(o;) r 1 (3.2) 

with self energy S(cj), is nondivergent (noncritical) even at T = below the mean field H C 2(0), 
because the dimensionality of the LLL fluctuation in 2D is two even at T = 0. Hence, it erases the 
mean field superconducting transition on the H C 2(T)-lme to create a first order vortex-solidification 
transition below it. In 3D, on the other hand, the dimensionality of dissipative LLL fluctuation 
at T = becomes three, and hence, a "ctitical field" H* 2 (0) can be defined at T = below the 
solidification field B m (0) (see §5 in ref.8). 

If the (Gaussian) VG transition occurs in a higher field than the disorder-free true transition 
(i.e., the solidification transition) point, a familiar resummation approximation of diagrams, such 
as RPA, is applicable© in adressing a glass transition point. Following previous works, we will 
invoke a systematic loop (or 1/M) expansion and focus on its lowest order (M = oo) terms by, 
as a mathematical tool, assuming ^> to have M-flavors. It is because this approach applied to the 
3 D thermal vortex state under a correlated (line-like) disorder parallel to the field has ledi to a 
transition line qualitatively consistent with the superconducting transition in high Tc superconduc- 
tors with parallel columnar defects. Note that, since the disorder with strength U p is, in eq.(3.1) 
at T = 0, persistently correlated along the "time" direction, the present quantum GL model with 
disorder is formally similar close to T = 0, to the bulk 3 D GL model at high temperatures with 
line-like defects parallel to the field. 

Hereafter, we write S(cj) = X vg (u;) + Sp, where the first and second terms are, respectively, given 
by Fig. 3 (a) and Fig. 3 (b). The expression of Xvg consistent with them has the form of a ladder of 
pinning lines with vertex correction, as described in Fig. 3 (c) and (d). In expressing the diagrams 
explicitly, the bare pinning strength U p / (2irr 2 B ) can be regarded as being replaced by 



= U P J^v k f 00 (k 2 r 2 B ) I 1 + a vg v k = ^C( C r vg ), (3.3) 



' B 

where the factor C(a vg ) = J °° alkk f 00 (k 2 )e~ k2 / 2 (l + a vg e~ k2 / 2 )- 2 with 



(7- 



vg 



E(^")) 2 ( 3 - 4 ) 



implies a renormalization (vertex correction) due to the LLL fluctuation of the pinning strength. 
The detailed form of C(o~ vg ) depends on the functional form of /oo an d thus, as mentioned in §2, 
is not known concretely. We just expect /oo(k) to have an algebraic form in k 2 . It implies that, 
when <7 vg S> 1, will roughly behave like ~ a~ 

s U p/ r B- The two self 

energy terms are given by 
S vg H = -A< ii) 6( W ) 1 (3.5) 
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* = |E/ k ^) = ^? eM - (3 - 6) 

and thus, eq.(3.2) yields 

( r 1 = /i(0) + 7 |w| - A^g{u;) + E F . (3.7) 

The Gaussian VG susceptibility Xvgi consistent with the above S vg ((j), is given as a series of ladder 
diagrams. Since the propagator Q(uj), due to the LLL degeneracy, depends only on a Matsubara 
frequency u, the sum of ladder diagrams is a trivial geometrical series, and hence, Xvg simply 
becomes 

Xvg = 27rr|A< i * ) ^l-A<V(0)j (3.8) 

so that the "Gaussian" VG transition point is given by the equation 

A^V(O) = ^f ( 2 Vg V (0) = 1. (3.9) 

Eq.(3.7) is quadratic in the renormalized propagator Q{oj) and easily analyzed. Further, at T = 
where the frequency summation is replaced by an integral, we can accomplish the present analytic 
calculation. Using eq.(3.9), the propagator Q(co) resulting from eq.(3.7) is expressed at the VG 
transition field -B V g,o by 



0-i( w) = 0-1(0) + + 



7H|a- 1 (0) + ™| (3-10) 



with 



2^ 1 (0)=ln(^ y l+S F , (3.11) 



where eq.(3.9) was used. The <7 vg -expression at -B V g,o is expressed as 



and, by combining this with eq.(3.9), a vg is a function of the combination 



and given as a solution of 



(3 - 13) 

<C(a vg ) = V . (3.14) 

On the other hand, the integral of requires a high-frequency cut off. By reasonably assuming a 
constant of order unity A c ~ jujm, where u>m is a frequency cut off, the frequency sum (integral) 
of eq.(3.6) results in 

-r = JtA^^ 1 <3 ' 15 ' 
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Applying this to eq.(3.11), we obtain the relation 

111 ^(0) J" ^r%,[3a vg ln [ U 4 ' (3-16) 

Eqs.(3.13), (3.14), and (3.16) give the T = VG transition field -B V g,o- I n the bracket of r.h.s. of 
eq.(3.16), the first term is a measure of the vortex pinning strength, while the second term is a 
measure of quantum LLL fluctuation effect. 

First, let us apply the above result to the ordinary dirty limit. According to eqs.(2.2) and (2.4), 
the coefficients in GL action take the following values in T —* limit 

7-7 (0) (T = 0) = ^, 
irB l z 

C/ 4 ^iV" 1 (0)(7 (0) (r = 0)) 2 , (3.17) 

together with eq.(2.10). Consequently, in this case, rj is a universal number, 5 x 10 -3 , and hence, 
<7 vg is also a constant, possibly of the order of 10 _1 . Further, the strength U4/ '(4-7r 2 r 2 ? 7) of quantum 
superconducting fluctuation becomes (2ttEft)~ 1 . Then, eq.(3.16) implies that, although -B V g,o will 
lie below H^ 2 (0) for large enough EpT, it increases with increasing the disorder strength l/{E-pr). 
One can find that the primary origin of this -B vg! o-increase is the disorder-induced enhancement 
of H^ 2 (0) oc t _1 , while the cancellation between the two terms in the bracket of eq.(3.16) is 
subtle and may depend on the diagram-resummation method which should be refined in higher 
orders of (Epr) _1 . Here we will respect the present result in the simple but consistent M _1 = 
approximation and expect, more generally (but in the dirty limit), the fluctuation contribution to 
outweigh the pinning contribution and to make -B vg ,o lower than H^ 2 {Q) at least up to 0((£ft) _1 ) 
(see also the next paragraph). On the other hand, without the fluctuation (second) term in eq.(3.16), 
it would show a T = superconducting transition point increasing with increasing disorder and 
existing above H^ 2 (0). This statement is essentially the same as the argument by Spivak and ZhouS 
of an unlimitedly large "i? C 2" (see also the sentence prior to eq.(2.7)). As already mentioned, we 
expect the inclusion of the quantum superconducting fluctuation to push the superconducting 
transition point down to -B vg ,o below H^ 2 {Q), even in the ordinary dirty limit. 

The true transition point B vg should be lowered further from the Gaussian one -B vg ,o by going 
beyond the present Gaussian approximation and taking account of interactions between the VG 
fluctuations. E X pl a iuiu g this suites auotfle, apparatus and wifl be g iveu in a separate paperB 
We note here that this shift 1 — B vg /B vg: Q is also of the order of (Epr) -1 . Nevertheless, this fact 
does not change the above statement in the ordinary dirty limit that B vg < H^} 2 (0), while B vg 
will increase with increasing (Epr) . Below, we will not distinguish B vg from -B V g,o m discussing 
effects of an electron-repulsion. This simplification does not affect conclusions which follow in this 
section. 
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Now, let us examine how an interplay between an electron-repulsion and disorder changes the 
above results in the ordinary dirty limit. According to the previous works© (see also ref.19), the 
contributions of the dynamically-screened Coulomb interaction to the linearized GL equation for 
quasi 2D films can be included perturbatively by assuming in 2D case a short-ranged repulsive 
interaction with strength Ai = R r j '(8ttRq) = 3/(kpdl), where d is the film thickness, and R r 
is identified with the high temperature sheet resistance. As shown in ref.9, the corrections due 
to the short-ranged electron repulsion to the GL coefficients /j,(0) and U4 are convergent in low 
T limit at each order in Ai, since the frequency dependences of Cooperons appearing as vertex 
corrections for the couplings to the pair-fields are cut off by the -B-dependence. It implies that 
/x(0) and U4 are determined by the high frequency side of Matsubara frequency summations and 
that their T dependences will, irrespective of Ai-values, be lost roughly when t < 1, i.e., T < 
T™ f = 0.15T c qB j H^ 2 (0) (see §2). As mentioned above eq.(2.7), the situation is also similar in 
Ai-dependences of C/ p /oo(krB). In contrast, each term of Ai-perturbation series for the time scale 7 
is logarithmically divergent, and, at low enough T, 7 has the systematic expansion parameter!© 
Ailn(T/T™ f ). It is expected from this systematic perturbation series that an onset temperature 
below which 7 begins to rapidly decrease on cooling will be given by 

T rep (Ai) ~ T£ { exp(-c rep /A 1 ), (3.18) 

where c rep is a positive constant (possibly) slightly less than unity. These facts imply that we 
have an intermediate temperature region below T™ f but above T rep (Ai). Besides these microscopic 
temperature scales, we have the quantum-thermal crossover temperature T cr = U4/ '(2-7rr^7 2 ) on 
the LLL fluctuation behavior (see §2 in ref.8). Then, if T rep < ^ (< T™ { ), there is a low but 
intermediate temperature region above T rep but below T cr in which any microscopic T-dependence 
carried by the GL coefficients is negligible and the LLL fluctuation is of quantum nature, although 
the GL coefficients are different from those in the dirty limit due to the Ai corrections. In §4, a 
computational evidence on the presence of such an intermediate temperature range will be given. 
Then, an apparent VG transition field B* g can be defined in this intermediate temperature region 
as a "T = 0" transition field and, according to eq.(3.16), is expressed by 

B* g ~H c2 (0)(l-d* s (E F T)~ l ), (3.19) 

up to 0((Sft") _1 ), where H C 2(0) and the constant d* g have Ai-corrections. Although we expect d* g 
to be positive, it is possibly smaller than unity, at least when Ai = 0, once recalling the above- 
mentioned cancellation between a fluctuation term and a pinning term. Although the Ai-correction 
to d g is due to those in U4 and U p , it is not easy to exactly obtain those corrections even up to 
O(Ai). However, since the Ai-correction in d* is accompanied in eq.(3.19) by the factor 1/EpT, 
the Ai-dependence of B* g can be seen as being dominated by that of H C 2(0). It is now known 
that H C 2(0) decreases!'© with increasing R r , and thus, the apparent critical field B* defined in 
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the intermediate temperature range, is expected to decrease with increasing R r , as indicated in 
experiments.BH!) 

To find the i? r -dependence of a true (but, possibly, inaccessible) B vg at T = 0, we need the 
Ai-dependences of the GL coefficients in low T limit T <C T rep . Since, at the present stage, we have 
no computation evidence enough to argue that 7 will remain positive in low T limit, we will merely 
assume here j(T — > 0) to approach a small but positive value 7 m i n . Then, according to eq.(3.16), 
B vg is given by eq.(3.19) with d* g there replaced by d g , which will take the form 

7 (0) / 7 (o) \ 

d g ^^—ln[?—\ »d* g , (3.20) 

6' 1 I mm \ I ram I 

where other 7 m i n -independent terms were neglected by assuming 7 m i n <C 7^ . Since it will be clear 
that 7min is insensitive to Ai or decreases with increasing Ai, B vg will decrease, more drastically 
than B* g , with increasing R r . However, we believe through the results in §4 and experimental 
informations that the field B c at which a flat (T-independent) resistance curve is seen should 
correspond to B* g and that B vg is not measurable because T rep is inaccessibly low. 

Finally, let us compare eq.(3.19) with the corresponding expression on B m (0). By substituting 
the parameter values in the dirty limit into the expressiorS |/xo(0)| oc C^4°V( 7rr B7^°^) derived in 
ref.8 (see §2.2 there), we obtain 

B m (0) ~ H c2 (0)(l - d^Epr)- 1 ), (3.21) 

where the constant d m was argued there& to be more than 6.0. By taking account of the prospects 
on <i*-value mentioned above, we believe here that, in general, B m (0) will lie below B* g . Further, 
due to the similarity on parameter dependences in eqs.(3.19) and (3.21), the difference between 
B m (0) and B* g should enlarge with increasing R r . 

§4. Computation of 7 at low T 

As mentioned in §3, we have previously^ shown that T(T) = 7(^/7^ (T = 0) at low enougn 
T takes a form of power serieslll' in Ailn(T/T c I J? f ) suggestive of the presence of a temperature 
scale T rep (Ai) of the form (3.18), where (T = 0) is the limiting value given in eq.(3.17). To 
demonstrate our argument in §3 on the presence of the T-insensitive intermediate region above T^gp, 
we have carried out a numerical calculation of T(T) useful even at low enough T on the basis of 
the resummation technique of Oreg and Finkel'stein (OF).© Its preliminary results will be briefly 
reported here. 

The OF's technique was originally developed to obtain the reduction of mean field T c in low 
dimensional s-wave case due to the interplay between the electron-repulsion and disorder and 
subsequently, was extended to nonzero field case in ref.19 to examine -ff C 2(T)-lines. This method 
focuses on the selfconsistent equation for the vertex part V c in the Cooper channel and consists 
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of solving it as a matrix equation in the (Matsubara) frequency space. The mean field transition 
point is obtained as a vanishing eigenvalue E c (u = 0) of the inverse of V c (uj = 0), where uj is 
the Matsubara frequency carried by the pair-field. One can manage to, in order to find T valid 
beyond the lowest order in Ai, extend this technique to the case with nonzero u> by applying^ a 
Pade-approximant to the frequency dependence of E c . 

In Fig4, the T v.s. T curves computed for various Ai values and at the field H C 2(0) of each Ai 
are given. In the temperature range T/T c q < 0.15 (i.e., T < T™ f ) of our interest, the temperature 
variations of T seem to become weaker with increasing Ai except at the lowest temperatures. The 
reduction of 7- values accompanying the Ai-increase at the intermediate temperatures arises, as well 
as that of i/ C 2-values, from the high frequency contribution in the diffusion propagators which was 
neglected for brevity in the 7-expression calculated in ref.9. It is not easy to judge how T rep should 
be defined from the curves, and it will be defined, for convenience, as the temperature below which T 
becomes less than T(T/T c q = 0.15). Then, we have, for instance, T Tep (R r / Rq = 0.257) ~ 0.005T c o, 
and T mp (R r /R Q = 0.428) ~ 0.013T c0 . Further, according to Fig.5 in which R r /R Q = 0.428 
is commonly used, the T rep thus defined decreases with decreasing B consistently with eq.(3.18) 
proportional to T™ f . Based on these figures, we conclude that the presence of an intermediate 
temperature region, in which 7 and hence, B* g can be defined as quantities insensitive to T, has 
been justified by the above microscopic computation. 

§5. Summary and Discussion 

In this paper, we have examined the position and parameter dependence of 2D VG transition field 
to be expected in low T limit. In the realistic model with repulsive interaction between electrons, 
we have a complicated situation due to a microscopic T-dependence of the time scale 7 becoming 
remarkable rather below a very low temperature scale T rep : An intermediate (but low) temperature 
range above T rep , rather than the low T limit below it, is expected to be relevant to experiments 
suggesting a 2D FSI behavior so that an apparent critical field B* g , being insensitive to T there, 
plays the role of a "T = 0" VG critical field. In a companion paper, we will show how available 
resistivity data suggesting a T = FSI transition can be explained based on the present results 
and argument. 

Finally, we wish to connect the present result with the resistive behavior near the disorder free 2D 
quantum melting linefl B m (T) ~ B m (0). In ref.8, it was pointed out that, at nonzero temperatures 
in the quantum regime T < T cr , the usual fluctuation conductance results in " fan-shaped" resistivity 
curves similar to the 2D FSI behavior near or below a fluctuation-corrected H C 2(0) and, only close 
to and below B m , reduces to the classical vortex flow behavior at the same T. Further, it was 
argued even that this itself may be the origin of the 2D FSI behavior. This argument based on 
the neglect of vortex pinnings may be justified only if B m > B* g . However, as seen at the end of 
§3, B m (0) is in general likely to lie below B* . Therefore, to try to understand comprehensively 
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the FSI behaviors in disordered thin films, a VG fluctuation contribution to conductance needs 
to be included which will be examined in ref.17. In fact, recent data analysis© has suggested a 
vortex lattice melting to occur much below the critical field B c at which the resistance becomes flat 
(see ref.17). We wish to note that, nevertheless, the knowledges!^ on the disorder- free fluctuation 
conductance at low T become important© in understanding differences in the intervening metallic 
resistance value at B* g between various materials. 

After submitting this manuscript, we were aware of the paper by Galitski and Larkinll) who have 
also examined a quantum transition field on a macroscopic superconductivity in disordered thin 
films from a different point of view and by neglecting an electron-electron repulsion. Although, 
in contrast to our approach, they have assumed a spontaneous creation of granular structure, 
nevertheless an expression determining the transition field at T = has been derived which is 
essentially the same as our eq.(3.16), except for the absence in their expression of eq.(3.16)'s second 
term (note that Eyr is denoted as g in ref.23), and has led them to a conclusion similar to our 
statement given below eq.(3.17) in relation to ref.14. However, the quantum amplitude fluctuation 
and the microscopic interplay between disorder and an electron-electron repulsion, both of which 
should contribute to a decrease of a quantum transition field, have been ignored in ref.23. Further, 
although a glass behavior is expected in lower fields at T = 0, the transition field was determined 
there, according to the sentences below their eq.(16), by assuming an occurrence of the ordinary 
phase coherence in contradiction to our theory with the amplitude fluctuation included showing that 
the glass ordering is not signaled by a development of the ordinary phase coherence. Explanation 
of experimental data listed in ref.23 should be ascribed to thermal fluctuation effects. 0© Details 
of this discussion will be given elsewhere.© 
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Figure Captions 



Fig.l 

Feynman diagrams contributing to T and belonging to the same family when classified according 
to how the diffusion propagators appear. Double solid lines imply the pair-field propagators, a single 
solid line is an electron propagator, the dotted line with open circle is a diffusion propagator between 
different replicas, the dotted line with cross denotes a single impurity line carrying (2-7riV(0)r) _1 , 
and the hatched corner vertex parts imply the Cooperons which modify the couplings to the pair- 
fields. 

Fig.2 

Details of the first three diagrams in Fig.l. 
Fig.3 

Diagrams necessary in obtaining B vg) o, where a hatched corner vertex is defined in (c) (which 
should not be confused with that in Fig.l and 2), the double dotted line denotes the pinning line 
carrying U p , and the solid circle is the interaction strength carrying U4. See the text for other 
details. 

Fig.4 

r(T) v.s. T curves computed in terms of OF's technique for different values of 87rAi = R t /Rq = 
(top), 0.086, 0.171, 0.257, 0.342, and 0.428 (bottom). Each curve was obtained by fixing B to the 
-ffc2(0)-value at each Ai-value, and 2ttT c qt = 0.25 was used commonly. 

Fig.5 

T(T) v.s. T curves obtained for different field values, B/H C 2(0) = 1.2 (top), 1.0, and 0.8 (bottom) 
by fixing 87rAi and 27tT c qt to 0.428 and 0.25, respectively. 
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